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ABSTRACT 

Aims. We have searched for temporal variations of narrow absorption lines in high resolution quasar spectra. A sample of five distant 
sources have been assembled, for which two spectra - either VLT/UVES or Keck/HIRES - taken several years apart are available. 
Methods. We first investigate under which conditions variations in absorption line profiles can be detected reliably from high reso¬ 
lution spectra, and discuss the implications of changes in terms of small-scale structure within the intervening gas or intrinsic origin. 
The targets selected allow us to investigate the time behavior of a broad variety of absorption line systems, sampling diverse environ¬ 
ments: the vicinity of active nuclei, galaxy halos, molecular-rich galaxy disks associated with damped Lya systems, as well as neutral 
gas within our own Galaxy. 

Results. Intervening absorption lines from Mg ii, Ee ii or proxy species with lines of lower opacity tracing the same kind of (mod¬ 
erately ionised) gas appear in general to be remarkably stable (Icr upper limits as low as 10% for some components on scales in 
the range 10 - 100 au), even for systems at Zab.i ~ Ze. Marginal variations are observed for Mgii lines toward PKS 1229-021 at 
Zabs = 0.83032; however, we detect no systems displaying changes as large as those reported in low resolution SDSS spectra. The lack 
of clear variations for low ji Mg ii systems does not support the existence of a specific population of absorbers made of swept up gas 
toward blazars. In neutral or diffuse molecular media, clear changes are seen for Galactic Nai lines toward PKS 1229-02 (decrease 
of A by a factor of four for one of the five components over 9.7 yr), corresponding to structure at a scale of about 35 au, in good 
agreement with known properties of the Galactic interstellar medium. Tentative variations are detected for H 2 7 = 3 lines toward 
EBQS J2340-0053 at Zabs - 2.05454 (== 35% change in column density, N, over 0.7 yr in the rest-frame), suggesting the existence 
of structure at the 10 au-scale for this warm molecular gas. A marginal change is also seen in C i from another velocity component of 
the same absorption system (=: 70% variation in A(Ci). 

Key words. Quasars: absorption lines - ISM: structure 


1. Introduction 

Absorption lines in distant active galactic nuclei (AGN) spectra 
have long been recognized as an invaluable source of informa¬ 
tion on diffuse material lying all along these lines of sight (LoS), 
from local interstellar material within the Milky Way, disks or 
halos of intervening galaxies, intergalactic clouds, up to gas lo¬ 
cated in the vicinity of the AGN itself. As illustrated by numer¬ 
ous studies (see e.g. Noterdaeme et al. 2007; Neeleman et al. 
2015) many properties of the gas can be investigated through 
a detailed analysis of the absorption lines, including ionization 
level, chemical composition, density and temperature. 

Aside from the broad absorption lines (BAL) and intrinsic 
narrow absorption lines which often display time variations over 
timescales of months or years in the absorber’s frame (Hamann 
et al. 2008, 2011; Filiz Ak et al. 2013 and references therein; 
Misawa et al. 2014; Grier et al. 2015) intervening absorption 
systems observed in quasar spectra are thought to be essentially 


* Based on observations with UVES on the Very Large Telescope at 
the European Southern Observatory (under programs 67.C-0157, 68.A- 
0170, 074.B-0358, 278.A-5048, 082.A-0569, 087.A-0597 and 087.A- 
0778) and with HIRES on the Keck Telescope. 


stable in time. However, a recent study based on multi-epoch 
low resolution SDSS spectra (Hacker et al. 2013) suggests that 
some narrow intervening systems (mainly Mg n and Fe ii) might 
also display variability. These authors propose that time changes 
are due to the transverse motion of the LoS through the absorber 
coupled to the presence of structure in the absorbing gas at scales 
in the range 10 - 100 au, similar to that observed for diffuse neu¬ 
tral gas in our own Galaxy (Crawford 2003; Lauroesch 2007, 
Welty 2007, Meyer et al. 2012). Indeed, transverse peculiar ve¬ 
locity values of a few 100 km s ' are expected for the target, the 
observer and the intervening gas, which imply drifts over a few 
years falling in this range (recall that 1 au yr * is equivalent to 
4.67 km s“'). 

Information on spatial structure within the low-ionization 
medium similar to the one probed by multi-epoch SDSS spec¬ 
tra remains quite limited. Observations of gravitationally lensed 
quasars with multiple images have been used to investigate 
spatial variations over scales of a few 100 pc. In particular, 
Rauch et al. (2002) studied the behavior of Mg 11 - Fe ii systems 
along three adjacent LoS toward Q2237I-0305 and by compar¬ 
ing the line profiles, infer a lower limit of about 0.5 kpc for 
the overall absorber extent and a typical size of the order of 
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100 pc for cloudlets associated with individual velocity compo¬ 
nents. Similarly, toward the lensed QSO pair APM 08279-1-5255, 
Kobayashi et al. (2002) derived an estimate for the cloudlet size 
of 200 pc, in good agreement with values inferred from mod¬ 
elling of these absorbers (Churchill & Charlton 1999; Crighton 
et al. 2015). Such observations are restricted to a few systems 
and do not enable a study of the structure at smaller scales. In our 
own Galaxy, structure in the low-ionization medium has been 
investigated by Welty (2007) who found no change in O i. Si ii, 
S II and Fe ii lines over tens of au. Thus, the limited information 
available (see also the recent study by Neeleman et al. 2015) 
does not suggest the presence of significant structure over scales 
as small as 100 au in the low-ionization medium, comparable to 
the one seen in the neutral medium. 

The results obtained by Hacker et al. (2013) then appear 
somewhat puzzling and require additional investigations. Their 
study is based on low resolution SDSS spectra which can only 
provide equivalent width measurements and do not give access 
to the true line profiles and to the identification of the individ¬ 
ual absorption velocity components. If some absorption lines 
do vary, changes most probably affect only some velocity com¬ 
ponents and should be much more apparent in high resolution 
spectra. High resolution line profiles would also be very use¬ 
ful to assess the reality of any variation detection. Indeed, when 
spectral profiles are resolved (or at least partially resolved), 
changes in a specific transition from a given species imply a cor¬ 
responding change for other transitions from the same species 
that can be computed a priori, providing a stringent reliability 
test. Furthermore, high resolution data will minimize the impact 
of various artefacts like imperfect sky subtraction, bad pixels etc. 

Good signal-to-noise (S/N) high resolution quasar spectra 
have been taken routinely over the last 15 years (for instance, in 
the context of Large Programs conducted on 8 - 10m class tele¬ 
scopes; see e. g. Bergeron et al. 2004 and Molaro et al. 2013) 
and for some distant sources, several observations have been 
performed over time intervals reaching ten years. We have thus 
gathered a sample of five sources with at least two VLT/UVES 
or Keck/HIRES spectra taken several years apart. In one case, 
we acquired a spectrum specifically for this study. These lines 
of sight altogether include a broad variety of systems, probing 
various environments ranging from neutral Galactic gas, galaxy 
disks at intermediate redshifts (i.e. damped Lya system). Mg ii- 
Fe ii systems due to intervening galaxy halos as well as systems 
with a redshift close to the emission redshift. Several of our tar¬ 
gets fall in the blazar category, which will allow us to investigate 
the possible presence of narrow systems with a large ejection 
velocity, in line with the scenario put forward by Bergeron et al. 
( 2011 ). 

The purpose of our study is to explore the potential of multi¬ 
epoch high-resolution observations of AGN to investigate the 
small-scale structure in intervening gas and possible intrinsic 
variations for narrow low ejection velocity systems. In Sect. 2, 
we describe the method (expected LoS drifts, selection of transi¬ 
tions which are best suited for a search of variations, etc) while 
the data and raw results are presented in Sect. 3. Their impli¬ 
cations are discussed in Sect. 4. We summarize our results and 
mention some prospects for further variability studies in Sect. 5. 

2. Detecting absorption iine variations 

When considering two spectra of the same target taken over 
some time interval, an important parameter is the drift of the LoS 
relative to a foreground absorber. We shall first discuss how one 
can estimate this scalelength over which the comparison of the 


two observations will provide structure information. Next, we 
examine which line parameters (the opacity in particular) yield 
the best sensitivity to a given fractional column density change. 
Finally, we discuss methods which can be used to reliably detect 
line variations and assess their significance. 

2.1. Line of sight drifts rotative to the absorber 

For now, we consider that the target is point-like and defines 
an “ideal” LoS, intersecting an intermediate redshift absorber 
at point M (hereafter the “impact point”; Fig.[TJ. As mentioned 
above, it is the combination of small-scale density or velocity 
structure in the intervening gas and the transverse motion of the 
LoS through this material which can potentially induce line vari¬ 
ations. The appropriate frame for computing the velocities of in¬ 
terest is a cosmological one, in which the Cosmic Microwave 
Background (CMB) is globally isotropic. The motion of both 
the target and the observer determines the drift of the LoS in this 
frame. The observer’s motion is due mainly to the peculiar mo¬ 
tion of the Milky Way and to Galactic rotation. The combination 
of the two corresponding velocities is directly constrained by the 
detection of the dipolar anisotropy in the CMB radiation indicat¬ 
ing that the Solar System is moving at 368 + 2 km s“' relative 
to the observable Universe toward the direction lo - 263.85 deg 
and bo - 48.25 deg (Fixsen et al. 1996). 

Let us denote by Vo the corresponding velocity vector. 
Consider a target T with Galactic coordinates {lT,bT) and let 
Ut be the unit vector pointing in this direction (Fig. [TJ. We are 
concerned only by the transverse component, thus by the projec¬ 
tion, Vo j., of Vo onto a plane normal to ux. 

This projected velocity can be written as 

Vo± = Vo - (ut - Vo)Ut = Vo (uo - (ux ■ uo) ut), (1) 

where Vo = Vo • Uo. From these relations, and using the compo¬ 
nents (cos bj cos Ij, cos bj sin Ij, sin bj) of the unit vector asso¬ 
ciated with Galactic coordinates [It, br), it is straightforward to 
compute Vo, j. for any specific target (for the five targets listed in 
Table 1, Vo,± = 288,368,129,365 and 12 km s“' respectively). 
In the most favorable case for which Vo,± = Vo, this motion 
alone results in a shift of about 4 x 10'^ pc or 800 au near the 
observer, over a time interval of 10 years. 

Unfortunately, we cannot estimate the transverse velocity of 
the target because the proper motion onto the sky is much too 
small to be measurable, contrary to interstellar medium stud¬ 
ies, in which the proper motion of nearby bright stars is known. 
Unless the target belongs to a peculiar environment like a galaxy 
cluster or experienced an ejection event, the magnitude of its pe¬ 
culiar velocity should be also of the order of a few 100 km s '. 
The transverse velocity of the impact point M, can then be writ¬ 
ten as, 

Vm^ =aVx^ + (l-a)Vo^, (2) 

where a = dAldj, d^ being the distance to the absorber and dj 
that to the target. Finally, one must take into account the velocity 
Va of the absorber in the cosmological frame or rather its trans¬ 
verse component, Vaj.. In the end, the value of interest is the 
relative drift velocity 

Vd = VM^-VAx. (3) 

Again, Va± cannot be determined (interstellar studies suffer 
from the same difficulty) and only a typical value of a few 
hundreds of km s * can be considered (as with the observer. 
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interstellar gas is at rest in average. In this case, the drift of the 
LoS is determined by the peculiar motion of the Sun in this 
frame. Vs, or more specifically, by its projection onto a plane 
perpendicular to the LoS considered. This question is discussed 
in more detail in Boisse et al. (2013). 


2.2. Sensitivity to absorption iine variations 

In order to select the absorption lines which are best suited to 
search for time variations, we now quantify the sensitivity to 
a given fractional column density (N) change. Obviously, very 
weak or, on the contrary, saturated features will be of little use: 
there must be some optimal intermediate value for the opacity. 

Consider a fully resolved normalised profile, /„(u), expressed 
in the velocity scale. 


Fig. 1. The observer and target move with 3D velocities Vo and 
Vt in the cosmological frame. Their transverse component Vo j. 
and Vtj. determines the transverse velocity Vmj. of the impact 
point M. The latter combines with the transverse velocity of the 
absorber Vaj. to define the relative drift velocity of the LoS with 
respect to the gas, Vd, which sets the scale probed within the 
absorber during a given time interval. 


the actual value results from the combination of the bulk inter¬ 
vening galaxy motion and of the absorbing gas internal kine¬ 
matics, e.g. rotation, infall or outflow). A maximum value for 
the drift velocity, Vj, can be estimated assuming an optimal 
configuration in which Vo± and Vjj. are parallel while Va± 
points in the opposite direction. In that case, Vd can reach values 
Vd = max{Vo±, Vt±) + Va±, as large as about 500 km s“'. 

Eqs. (2) and (3) are valid in the local 3D space. On large 
cosmological scales, there is no unique way to define distances. 
Thus, these relations can hold only at low redshifts where all 
cosmological distances become equivalent. Since the exact val¬ 
ues of both Vt ± and Va ± are unknown for any specific target, it 
is sufficient for our purpose to estimate drift velocities by i) in¬ 
serting into the above equations typical values for Vt± and Va±, 
ii) using the angular size distance for d, and iii) taking into ac¬ 
count time dilation (a time interval At in the observer’s frame 
corresponds to Af/(1 +za) or Af/(1 -i-zt’) in the absorber’s or tar¬ 
get’s frame respectively, which is equivalent to reduce Va± and 
Vt-j. by a factor of (1 H- za) or (1 -H zr))- 

In fact, a quasar is an extended continuum source with a typ¬ 
ical size in the range 10 - 100 au (the size of the accretion disk 
responsible for the optical emission; Dai et al. 2010). Thus, con¬ 
trary to observations of moving stars used to probe interstellar 
gas structure, significant averaging of the foreground structure 
can occur, implying a reduced amplitude for absorption line vari¬ 
ations. The situation may be different for blazars. These sources 
undergo spectacular photometric variations and, when in a “high 
state”, their continuum is dominated by emission from the rela¬ 
tivistic jet. Changes in the morphology of the source itself can 
also occur, which would be equivalent to an additional contribu¬ 
tion to the drift of the LoS. Time changes in the 21 cm absorp¬ 
tion seen toward AO 0235 h- 164 and PKS 1127-14 by Wolfe et al 
(1982) and Kanekar & Chengalur (2001) are precisely attributed 
to such an effect. 

When dealing with absorption by Galactic interstellar mate¬ 
rial (four out of our five selected targets display relatively strong 
Can or Nat absorption from Milky Way gas), the appropriate 
frame to consider is the Local Standard of Rest (LSR) in which 


IM = (4) 

where t{u) is the line opacity at velocity v, written as 


T{v)^kAofN(v), (5) 

where k is a constant. To the rest wavelength of the transition, / 
the oscillator strength and N{v) dv the column density for species 
within the velocity interval [v, v + dv]. Let us first assume a spec¬ 
trum noise which is independent on /„. The sensitivity can then 
be expressed as 


5In 

6N/N 


- T e 


(6) 


It reaches its maximum at t = 1 for which /„ 0.37. The interval 
over which the sensitivity is reduced by less than 20% is 0.61 < 
T < 1.53, which corresponds to 0.22 < /„ < 0.54. 

An alternative case of interest is that of photon noise, for 
which the spectrum r.m.s. follows a ^fT^ law; the ability to detect 
a given variation Sin is then proportional to 


din 

yfhrSNIN 


= re 


(7) 


This expression reaches its maximum at r = 2 where /„ ~ 0.14 
(the sensitivity remains good in the sense defined above for 
1.2 < T < 3.1). In practise, observed spectra are within these 
two extreme cases; we conclude that the sensitivity for the de¬ 
tection of variations is optimal for intermediate opacity values 
(typically 1.5), corresponding to normalised intensities of about 
0 . 2 . 

One should keep in mind that generally, depending on the 
instrumental resolution and Doppler b parameter value, absorp¬ 
tion profiles will not be fully resolved, implying apparent opaci¬ 
ties lower than the true ones (Savage & Sembach 1991). For in¬ 
stance, with an instrumental resolution of 6.6 km s ', the depth 
of a Gaussian line is reduced by a factor of 0.45 for b - 2 km 
s '. Then, if the true opacity at line center is tq = 1, the apparent 
opacity becomes 0.33. Therefore, to conclude on the sensitiv¬ 
ity of a given line to changes in N, one cannot rely only on its 
observed profile. Fortunately, several transitions from the same 
species are often available with various / values. Simultaneous 
fitting of the corresponding line profiles will then allow us to 
determine the true opacity and assess whether or not these tran¬ 
sitions are appropriate for a variability search. 
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2.3. Characterisation of line variations 

Our aim is to detect line profile variations and to quantify their 
significance. A first direct indication can be obtained by super¬ 
posing the two successive spectra after careful normalisation and 
adjustment of the wavelength scales if necessary. Such a raw 
comparison is meaningful only if the spectral resolution, R, is 
the same at the two epochs. If this is not the case, the higher 
R spectrum has to be degraded to the lower resolution value. 
Fortunately, the HIRES and UVES data used in this paper have 
very similar nominal resolutions - FWHM = 6.25 and 6.6 km s ' 
respectively - and we could verify that, given the S/N ratio, such 
a slight difference has a negligible effect on the profile of the ab¬ 
sorption features discussed in this paper (cf Sect. 3). Differences 
in the line spread function (LSE) might also play a role. One can 
use unresolved lines to check that such an effect does not affect 
the comparison of successive profiles. 

It turns out that for the systems investigated here, the two 
successive spectra available rarely display marked, unambigu¬ 
ously significant differences. Thus, in order to i) increase our 
ability to detect weak variations by relying on the presence of 
several transitions from a given species (e.g. Eeii, or H 2 ) and ii) 
take into account differing R values if necessary, a more pow¬ 
erful method is to simultaneously fit all profiles corresponding 
to these transitions, at each epoch. To this purpose, we used 
vpfitE] a routine which uses multicomponent Voigt profiles 
(convolved with the instrument profile, which we assume to be 
Gaussian here) and yields the redshift z,, column density N, and 
Doppler parameter for each velocity component, together with 
their Icr uncertainties. Comparison of the values obtained for 
each epoch then allows us to determine whether some of these 
parameters have undergone significant time changes. In all cases, 
the same sets of components were appropriate to fit the profiles 
at both epochs, with no significant change in Zi, nor in f>;. Some 
time variations in V, are detected or suspected for a few absorp¬ 
tions systems (see Sect. 3). To assess their significance, we com¬ 
pute the difference between epoch 1 and 2, AN, = V ,,2 - 
and its uncertainty cr(AV,) by combining uncertainties on V,;i 
and V, 2 in quadrature. This procedure provides a formal signif¬ 
icance which properly takes into account the noise level in both 
spectra. Other sources of uncertainty related to continuum place¬ 
ment or to the choice of velocity components are more difficult to 
quantify, especially when variations are small. As a conservative 
approach, variations with a significance level above 3.5cr and an 
associated relative column density variation, AN/{N) above 25% 
are qualified as ’’significant”. Variations at more than 2cr but 
which do not fulfill these two criteria are considered as marginal. 

When several transitions from the same species are avail¬ 
able, the quality of the fit obtained provides a stringent test of 
the internal consistency of the data and gives a powerful means 
to rule out false variations due to instrumental artefacts. Indeed, 
changes must not only be present for all transitions but their 
magnitude must be consistent with the / values. In this regard, 
species such as Feii, Ci, or H 2 , which display many transitions 
from a given level within a relatively narrow wavelength inter¬ 
val are especially suitable to probe weakly ionized, neutral or 
molecular gas, respectively. 

Often, several species can be used to trace a given phase 
(e.g. Mg II, Ee ii, S ii for weakly ionized gas). One can then select 
among all lines detected from these species those i) which opac¬ 
ity falls in the optimal range discussed above, and ii) occurring 
in spectral regions of good S/N ratio. 


' http://www.ast.cam.ac.uk/ rfc/vpfit.html 


3. Variability study of five blazars with 
low-moderate p systems 

3.1. The absorber sample 

One of our initial motivations was to explore the link between 
gas swept by the AGN jet and low (3 absorbers, as suggested 
by Bergeron et al. (2011). We then investigated among their 
blazar sample objects with jets, thus radio-loud and y-ray emit¬ 
ters (detected by Egret, Integral or Eermi/LAT). An additional 
selection criterion is the existence of Mg ii absorption systems 
at low relative velocity y6 = Av{Zsrn,Zab.T)lc < 0.15 with a rest- 
frame equivalent width W, (2796) > 0.3 A, as well as intervening 
systems and Damped Ly-a (DLA) systems at high fi. We then 
searched the VLT-UVES and Keck-HIRES archives for high- 
spectral resolution spectra of this blazar subsample, and selected 
targets observed several years apart, possibly about one decade. 
We checked that the signal-to-noise per pixel was S/N >10 for 
the main absorptions of interest: the Mgii doublets at y6 < 0.15 
and the Galactic Nai doublets. 

One blazar satisfying the above selection criteria, is the ex¬ 
tensively studied, highly variable BE Lac object AO 0235-1-164 
for which variations of 21 cm intervening absorption at Zabs = 
0.524 had been reported three decades ago (Wolfe et al. 1982). 
This absorber has been identified as a BAL quasar with very 
faint extended emission (Burbidge et al. 1996). The variations 
were interpreted as most probably due to variations in the bright¬ 
ness distribution of the background radio source. Another blazar 
is PKS 1741-038, which was re-observed by us with UVES in 
2011 to search for variability of its two strong Mgii absorbers at 
P - 0.074 and 0.287. Eor both objects, the available UVES and 
HIRES spectra have been obtained at epochs about 9.5 yr apart 
in the observer’s frame. 

We also searched for high resolution spectra of PKS 
1229-02, a blazar for which we obtained deep HST imaging 
showing optical emission associated with knots in the radio jet 
(Le Brun et al. 1997). There is a DLA system at Zabs = 0.3950 
with associated 21 cm absorption, and several Mgii absorbers at 
>6 < 0.15. The above selection criteria are also met for this source 
and there are two sets of high S/N UVES spectra taken at epochs 
9.2 yr apart. 

Einally, there are HIRES and UVES data available for two 
targets with Ci absorption at Zabs ~ 1-5 with associated 21 
cm absorption at high p (Kanekar et al. 2010): PKS 0458-020 
(^^bs = 1.5606) and FBQS J2340-0053 (Zabs = 1.3608), the lat¬ 
ter with observations about two years apart. In the radio spec¬ 
trum of PKS 0458-020, variable 21 cm absorption on time- 
scale of months has been detected at Zabs = 2.0395 (Kanekar 
et al. 2014 and references therein). There are also DLA systems 
at low p < 0.10 toward PKS 0458-020 at Zabs = 2.0395 and 
FBQS J2340-0053 at Zabs = 2.0545. The physical conditions 
in the low p DLA system of EBQS J2340-0053, with associated 
CI absorption, have been investigated by Jorgenson et al. (2010). 
The multiple system at Zabs = 1.3608 toward EBQS J2340-0053 
has been extensively studied by Rahmani et al. (2012) to con¬ 
strain the variation of fundamental constants. 

Information on the sample and the observations is given in 
Table 1 and that on the absorption systems in Table 2. Eor some 
UVES data, we used the reduced data sets produced internally 
at ESO, based on the latest calibration pipelines and the best 
available calibration data. Eor the other UVES observations, the 
spectra were reduced by us or made available to us by several 
colleagues (see the acknowledgements). Eor HIRES spectra, the 
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Table 1. The sample. 


target 

common name 

class 

coordinates 

J2000 

^em 

spectrograph 

date 

A range 
nm 

At, n 
min 

AO 0235+164" 

BL-7 

023838.9+163659 

0.937 

HIRES 

12-1997 

367-612 

40,4 





UVES 

02-2007 

303-388+58-668 

60,4 

PKS 0458-020 

FS-r 

050112.8-015914 

2.286 

HIRES 

02-1995 

394-629 

480" 





UVES 

10-2004 

390-665+670-850 

75,4 

PKS 1229-02 

FS-r 

123200.0-022404 

1.044 

UVES 

02-2002 

326-445+458-668 

60,3 





UVES 

04-2011 

326-445+458-668 

50,1 

PKS 1741-038 

FS-r 

174358.8-035004 

1.054 

UVES 

06-2001 

414-621 

40,2 





UVES 

04-2011 

414-621 

50,5 

EBQS 12340-0053 

FS 

234023.7-005327 

2.085 

HIRES 

08-2006 

306-589 

250" 





UVES 

10-2008 

330-450+478-681 

75,6 


class - BL: BL Lac object, FS; Flat Spectrum Radio QSO (FSRQ), y: Gamma-ray emitter. The four y-ray emitters are blazars. 
spectrograph - spectral resolution: FWHM(FtIRES)=6.25 km s“', FWFtM(UVES)=6.6 km s“*. 

Af, n Exposure time per individual exposure and number of exposures. 

^ In a very high state at both epochs. The new is derived from [Ne iv] emission. 

“ Total exposure time. 


data sets were reduced by Prochaska’s team. All redshifts and 
velocities considered in the following are vacuum heliocentric. 

In the analysis presented below, we estimate the variation of 
N between the two epochs and use the combined quadratic errors 
to evaluate its significance level, as described in Sect. 2.3. In case 
of no variation, we give the Icr limit on AN/{N), restricting this 
evaluation to fairly good S/N, thus meaningful limits: AN/{N) < 
50%. 

3.2. Tracers of the Interstellar Medium 

The main tracers of the local Interstellar Medium (ISM) de¬ 
tectable in our spectra are the Nai and Can doublets, Ca 1/14227 
as well as the CH^/14232 and CH/14300 molecular lines. Diffuse 
interstellar bands (DIBs), mainly at 4428, 5780 and 5797 A, 
could also be identified in some UVES spectra. 

Toward AO 0235-1-164, the Galactic Nat absorption shows 
a strong, single component at Zabs = H-0.00002. A simultaneous 
fit to both members of the Nat doublet yields identical values 
at the two epochs, 9.2 yr apart, within the Icr errors; A(Nai) 
= (2.51+0.04) xlO'^cm-^. 

The Galactic Nai absorption toward PKS 1229-02 is mul¬ 
tiple (5 components) and extends over 76 km s '. Significant 
variation is present for the component at Zabs = -0.000127: in 
9.7 yr, it has nearly disappeared, with A(Nai) decreasing by a 
factor of four (see Table 3 and Fig. |2). The other three strong 
components (A(Nai) > lO" cm“^) did not vary (Icr; 17%). 

Toward PKS 1741-038, Galactic molecular absorptions at 
wlsr = 2 km s“' have been extensively investigated, with de¬ 
tections of HCO^ (Lucas & Liszt 1996), CO (Akeson & Blitz 
1999) and C 2 H (Lucas & Liszt 2000), together with variable 21 
cm absorption (Lazio et al. 2001). The associated Na i absorption 
at Zabs = -0.00033 is multiple (4 components), partly saturated 
and blended with atmospheric Nai emission, thus preventing a 
search for variations. Two well-known DIBs are detected at both 
epochs, at 5780 and 5797 A. Their equivalent widths are equal to 
0.33 and 0.10 A respectively, and are similar at two epochs 9.8 
yr apart, with AWI{W} =21% +13% and 20% +23%. Two other 
DIBs at 4428 and 6196 A, the Can doublet and CH^/14232 are 
detected but are far too noisy in the 2001 spectrum to investigate 
variability. 


Finally, toward FBQS J2340-0053, the Nai 5897 line is 
not covered in the HIRES spectrum. The single component of 
the Can doublet at Zabs = -0.00002 remains stable (Icr; 3.8%) 
within 1.9 yr. 


3.3. Neutral and molecular gas at high z 

Toward PKS 0458-020, Ci is well detected at Zabs = 1.5605 
in the HIRES spectrum (with a single component, as reported in 
Kanekar et al. 2010) but of too low S/N in the UVES spectrum to 
get useful variability constraints. In FBQS J2340-0053 spectra, 
Ci is detected at Zabs = 1.3606 and 2.0545. The three well sep¬ 
arated components of the lower z system are stable (Icr; 12%). 
The fit of CI at Zabs = 2.0545 involves a minimum of eight com¬ 
ponents, two sets being heavily blended (components 3-4 and 
5-6 respectively). Within the blend at Zabs = 2.05473 a very nar¬ 
row component (number 6) shows a possible variation AN/{N) 
of 69% in 0.71 yr (rest-frame) at the 3.0cr significance level (see 
Table 3 and Fig. [3] left panel); the Icr limits for the other seven 
components are in the range 9-50%. 

Molecular hydrogen has already been detected toward 
FBQS J2340-0053 at Zabs = 2.0545 (Jorgenson et al. 2010). 
The Lyman and Werner transitions display six components 
(which roughly correspond to C i components number 1, 2, 3-4, 
5-6, 7 and 8) and there is a possible variation of the sixth H 2 , 
7 = 3 weak component at Zabs = 2.05513. Using the unblended 
components of the transitions at 1081, 1084 and 1099 A we 
find a decrease, AN/{N} = -34%, in only 0.71 yr (rest-frame), 
significant at the 3.Icr level (Table 3 and Fig. [3] right panel). 
Unfortunately, the corresponding Ci features are too weak 
to check whether Ci lines display similar variations (note in 
the left panel of Fig. |3]that features at t; = 40 km s ' for the 
C1/11277 and C1/11328 transitions are in fact dominated by 
CT transitions from the strongest components). Conversely, 
the fourth H 2 component which corresponds to the second C i 
blend (components number 5 and 6) with marginal variations is 
opaque in the available transitions and thus, insensitive to time 
changes. 
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Table 2. The absorption systems. 


target 

^abs 

P 

main 

doublet 

strength 

S/N“ 

associated low ions 
present in the 2 spectra 

AO 0235+164 

+0.00002 


Nai 

w 

67-56 

none 


0.52432 

0.235 

Mg II 

VSm 

28 - 14'’ 

Mg I, Fe II, Mn ii, Ca ii, Ti ii 


0.85142 

0.045 

Mg II 

w 

62-50 

none 


0.85237 

0.045 

Mg II 

s 

62-50 

Mg I, Fe II 


0.85562 

0.043 

Mg II 

W 

62-50 

none 

PKS 0458-020 

0.89025 

0.503 

Mg II 

s 

11 - 13 

Fell 


1.52712 

0.257 

Mg II 

VSm 

15*-5 

Fell 


1.56055 

0.244 

Mg II 

s 

14'’ - 5 

C I, Fe II, Zn ii, Cr ii, Ni ii 


2.03956 

0.078 

Mg II 

s 

10''' - 5 

Fell, Znii, Cm, Niii, Sin 

PKS 1229-02 

-0.00016 


Nai 

W 

52- 15 

Can 


0.39458 

0.364 

Mg II 

VSm 

44-13 

Mg I, Fe II, Mn ii, Ca ii, Ti ii 


0.69990 

0.182 

Mg II 

Vw 

50-13 

none 


0.75643 

0.150 

Mg II 

WS 

50-13 

Mg I, Fe II 


0.76836 

0.144 

Mg II 

w 

50-14 

none 


0.83032 

0.110 

Mg II 

w 

53-14 

Fen 


0.83102 

0.110 

Mg II 

Vw 

53-14 

none 


0.83128 

0.109 

Mg II 

w 

53-14 

Fen 

PKS 1741-038 

-0.00070 


Nai 

w 

12-32 

none 


-0.00033 


Nai 

S 

12-32 

DIBs 


0.52739 

0.288 

Mg II 

w 

5-13 

none 


0.52822 

0.287 

Mg II 

S 

5-13 

Mg I, Fe II 


0.67846 

0.199 

Mg II 

w 

8 - 18 

none 


0.90799 

0.074 

Mg II 

VSm 

10-28 

Mg I, Fe II 

FBQS J2340-0053 

-0.00002 


Nai 

S 

29-68 

Can 


1.36064 

0.261 

Mg II 

VSm 

30* - 64 

C I, Fe II, Zn n, Cr ii. Si n, Ni ii 


2.05454 

0.010 

Fell 

VSm 

34*' - 56*' 

Ci, Ci*,Feii, Sn, Sin 


“ : S/N per pixel around the main absorption for the older and recent epoch, respectively. 

: Mg II region not covered; we list instead the S/N around Feii2374-2382, Feiil608-1611. 
VSm : very strong, highly multiple absorption, »;r(2796) > 2.0 A. 

S : strong absorption, w^{2196 or 5891) > 0.3 A. 
w ; weak absorption, 0.02 < Wr(2796 or 5891) < 0.3A. 
ws : although weak, the Mg ii doublet is saturated. 

Vw : very weak absorption, 10^(2796) < 0.02 A. 


3.4. The Mgii-Feii systems 
3.4.1. Low (3 Mg it systems 

There are only two blazars with low /? Mg ii systems of moder¬ 
ate strength. For the other three, these systems are either DLAs 
or extremely strong saturated Mg it absorbers; we have then 
searched for proxies of intermediate opacity (e.g. Feii lines) to 
constrain variability. 

Toward AO 0235-1-164, there are three Mgii systems at low 
P'. Zabs = 0.8514 (three fitted components), 0.8523 (seven com¬ 
ponents) and 0.85562 (single). There is no variation of the triple 
and single systems (Icr limits on ANI{N} of 9.5% and 6.3%, 
respectively). For the multiple system, one Mg ii component at 
Zabs = 0.85236, blended with two weaker ones at Zabs = 0.85231 
and 0.85241, shows a tentative variability of its column den¬ 
sity AN/{N) of 53% in 5.0 yr (absorber rest-frame), but at a 
low significance level of only 2.3cr (see Table 3 and Fig. |4]i. 
For these three components, the associated Fe ii doublet is sta¬ 
ble (Icr: 5.2%) 

Toward PKS 1229-02, there is a triple, unblended system at 
Zabs = 0.830. The strongest component at Zabs = 0.83032 shows 
a marginal small variation AN/{N) of 17% in 5.0 yr (rest-frame) 
at a 3.1cr significance level (Table 3 and Fig.|5ll. There is no vari¬ 
ation of the associated four strongest Feii lines (Icr: 7.3%). The 
systems at Zabs = 0.83102 and 0.83128 are too weak to derive 


meaningful limits on their variabilty. There are two other Mg ii 
systems atyS < 0.15, and they are stable: the weak, single compo¬ 
nent doublet at Zabs = 0.7684 (Icr: 8.7%) and the triple system at 
Zabs = 0.7564 with one nearly saturated component (Icr: 20%). 

Only one of the three targets with extremely strong saturated, 
low p Mg II systems has proxy lines yielding meaningful lim¬ 
its. Toward PKS 0458-020, the Niii triplet and the Znii-Crii 
lines at Zabs = 2.03956 are all of moderate strength and heav¬ 
ily blended without clear variation. In PKS 1741-038 spectra, 
the Fell lines are saturated. Thus, some variability constraint 
can only be placed on the Zabs = 2.0545 DLA absorber to¬ 
ward FBQS J2340-0053: the three strongest S ii components of 
medium strength at Zabs = 2.0546 and the strongest one of Aim 
at Zabs = 2.05511 do not show variability (Icr: 32% and 28%, 
respectively; see Sect. 4.2). 

3.4.2. High p Mg ii systems 

For the high p systems, variations are not detected from either 
their Mgii doublets or proxy lines. Toward AO 0235+164, the 
two strongest components of the Can doublet at Zabs = 0.5239, 
associated with the DLA, do not vary (Icr: 26%). Toward PKS 
0458-020, the Feii/l/l2344,2382 lines at Zabs = 1.5271 (select¬ 
ing the 3 unblended components) do not vary either (Icr: 15%). 
Finally, toward FBQS J2340-0053, the highly multiple Fe ii sys- 











Velocity (km/s) Velocity (km/s) 


Fig. 2. Galactic Nai absorption toward PKS 1229-02: spectrum, its error (red curves) and simultaneous fit to both transitions (blue 
curve). The UVES 2002 data are shown in the left panel; in the right panel we display the UVES 2011 data together with the fits to 
the 2002 data (blue). Note the disappearance of the component at iiheiio = 9.6 km s“' (thick vertical tick mark), within a period of 
9.2 yr. The heliocentric redshift z = -0.000159 (strongest component) corresponds to r = 0 km s '. 


tern at Zabs = 1.3606, associated with Ci and 21 cm absorp¬ 
tions, is stable (Icr: 14% for all the components, using the weak 
Ee 11/1/12249,2260 doublet to constrain the other nearly saturated 
lines). 

4. Discussion 

4. 1. Presence of intrinsic systems in blazar spectra ? 

Time variability of narrow absorption lines at low yS have been 
securely detected in small samples of high ionization systems 
(involving C iv, N v and O vi), observed at low and intermedi¬ 
ate resolution (Narayanan et al. 2004; Wise et al. 2004, Hamann 
et al. 2011). The lines showing variations are essentially from 
mini-BALs and associated {Au < 5000 km s ') narrow-line 
systems. The time-scales investiged (in the quasar rest-frame) 
range from 0.28 to 8.8 yr, and the variable narrow lines remained 
fixed in velocity. High and low ionization systems (at fi > 0.01) 
with repeated observations have been identified by Hacker et 
al. (2013) in a very large sample of SDSS-DR7 quasar spec¬ 
tra: 33 out of 1084 systems have absorption lines with variations 
in their equivalent widths. Eor most of them Wr > 0.3 A and 
AWrjWr > 50%, with two-thirds of the variable systems being at 
’relatively low’ p (< 0.22). Low ionization species are present in 
all the 33 variable systems and C iv is detected in the 17 systems 
at sufficiently high Zabs (> 1-50). Chen et al. (2013) focussed 
their variability study on associated {fi < 0.03) Mgii systems 
identified in SDSS-DR7 and SDSS-DR9 (BOSS survey) spectra. 
They do not detect time variations at a significance level larger 
than 3cr in their sample of 36 Mg ii doublets. 

No variation of intervening Mgii systems in high spectral 
resolution spectra has been reported so far. Among the five tar¬ 
gets in our sample, three blazars have in total six Mgii dou¬ 
blets at < 0.15 which are suitable (unsaturated, not heav¬ 
ily blended lines) for variability studies. As described in Sect. 
3.4, only marginal column density variations have been detected 


for two Zabs ~ 0.85 Mg II absorbers at a significance level of 
2.3cr and 3.7cr, respectively. The other four Mgii absorbers, at 
2abs ~ 0.75-0.85, are stable with ANI{N) in the range 6-20% at 
Icr confidence level. Among them, the one in the spectrum of 
AO 0235-b 164 is 150 km s“' away from the tentatively variable 
system. 

In Bergeron et al. (2011), we found a potential excess of 
Mg II systems at /3 ^ 0.1 in blazar spectra and proposed that 
gas swept up by the AGN jet could be responsible for it. In such 
a scenario, absorption line variability can be expected, but, in 
the absence of predictions from numerical simulations, it is dif¬ 
ficult to anticipate the time-scale and magnitude expected for the 
column density changes (this should strongly depends on the oc¬ 
currence of marked variations in the activity of the nucleus itself 
at about the time when the blazar was observed). Thus, we can¬ 
not really draw conclusions from the lack of clear evidence for 
time variations among the limited set of systems which we could 
investigate here. 

4.2. Structure in moderately ionized halo gas 

We now discuss systems involving narrow absorption lines asso¬ 
ciated with halo gas which is optically thick to UV ionizing radi¬ 
ation, i.e. gas classically traced by Mgii or Eeii transitions. We 
mentioned above the few cases for which we could get stringent 
upper limits on column density variations for discrete velocity 
components seen through intermediate opacity lines. One good 
example is the Zabs = 1.3606 system toward EBQS J2340-0053 
involving 10 Eeii components. Even when most of the absorp¬ 
tion profile is optically thick (this is often the case for the Mg ii 
doublet), the data can be used to search for variations near the 
profile edges (e.g. appearance or disappearance of faint compo¬ 
nents). We could not identify any change of this kind. Moreover, 
since the component redshifts are let free when fitting the pro¬ 
files, we could also get constraints on velocity variations; again. 
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Fig. 3. C 1^1277,1280,1328 and molecular Hj, 7=3 71081,1084,1081 absorptions at = 2.05472 toward FBQS J2340-0053 in 
the November 2006 (HIRES, black line) and October 2008 (UVES, red line) spectra. Left panel: the eight C i components required to 
fit the profiles are shown by blue tick marks and the marginally varying feature is at u = 0 (thick tick mark). In the top C i71280 plot, 
three Ci71656 lines from the Zats - 1-3606 system are also shown (green tick marks). For the Ci71277 and Ci71328 transitions, 
we show the position expected for C i* transitions associated with the 8 components (cyan tick marks); note that the features at u ^ 
40 km s * are mainly due to C i* transitions. Right panel: six H 2 7=3 components are required to fit the profiles. Marginal variability 
is seen in the sixth component (thick tick mark) at u = 40 km s ' (some difference is also seen in the fifth component but it is 
significant only at the 2.0cr level, instead of 3.1cr for the sixth one). For both narrow components displaying marginal variations in 
Ci (u = 0) or H 2 {v ^ 40 km s '), we checked that the slightly higher resolution of the HIRES spectrum cannot account for the 
observed difference between the profiles. 


we found no significant change, with 3cr upper limits generally 
below 1 km s '. Occasionally, some line profiles are observed 
to be rather smooth as seen for the resolved S ii71253 and 1259 
lines at Zabs = 2.05454 in the FBQS J2340-0053 spectrum. As 
shown in Fig. |6] no significant difference is seen along the en¬ 
tire profile. The Icr upper limit on the relative flux variation is 
1.7 10“^ per resolution element along the S n71259 line profile. 
Since the latter spans over 120 km s ', the previous upper limit 
applies to about 120/6.6 = 18 distinct resolution elements. 

To conclude for these systems, we detected no structure in 
the density and velocity field within Mg ii/Fe ii halos at the scales 
that we could probe (typically of the order of 100 au). This is 
consistent with previous observations of moderately ionised gas 
in galaxy halos. Indeed, both Rauch et al. (2002) and Kobayashi 
et al. (2002) find that a typical value for the cloud size associated 
with the individual velocity components is 200 pc, implying very 
little structure at the much smaller scales that we probe here. 
We also note that, from the analysis of the number density of 
strong Mg ii systems detected on top of quasar emission lines in 
SDSS spectra, Lawther et al (2012) derive a lower limit of about 
0.03 pc for the individual Mg ii absorbing clouds. Finally, Welty 
(2007) analysed multi-epoch HST spectra of HD 219188 and 
detect no variations for S 11 , O i, Si ii and Fe ii in the disk of our 
own Galaxy, over scales of the order of 100 au. Thus, our results 
are in good agreement with all known constraints on structure in 
moderately ionised gas. 


In this context, the variations tentatively detected by Hacker 
et al. (2013) for narrow intervening Mg ii systems are most sur¬ 
prising. These authors analysed more than 1000 reliable quasar 
absorption systems with multi-epoch observations and extracted 
a small subset of 33 systems displaying some evidence for time 
changes in their equivalent widths. Their study differs markedly 
from ours in the number of systems investigated and, given the 
small fraction of systems showing variability - about 3% - it 
is true that we expect very few variable lines in our sample. 
However, it should be stressed that the two studies differ also 
strongly in the sensitivity expressed in terms of AWr or AWrlWr- 
The variations reported in the SDSS spectra are very large with 
in most cases AVI/. > 0.3 A and AVV/-/VV/. > 50%, correspond¬ 
ing to even larger fractional column density variations. For com¬ 
parison, the marginal variation quoted above in the Mg ii72796 
line at Zabs = 0.83032 toward PKS 1229-02 is smaller than 
AW, = 0.01 A and AW,/W, = 5%. 

If the large variations quoted by Hacker et al. (2013) are real, 
one can imagine that the number of systems expected to display 
changes at the level that we can reach in our high resolution 
study should be much larger, coming in contradiction with the 
absence of clear detection of variations in our sample. Before 
drawing any conclusion, additional observations aiming at con¬ 
firming the status of the 33 potentially variable systems would 
be highly desirable. In particular, since variations are (to first or- 
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Fig. 5. Intervening Mgii absorption at Zabs = 0.83032 toward PKS 1229-02: spectrum and its error (red curves) together with 
simultaneous fit to both transitions (blue curve). In the left panel, we show the UVES 2002 data, while the right panel displays the 
UVES 2011 data with the fit performed on the UVES 2002 spectrum. The blue vertical tick mark underlines the marginally variable 
Mg n doublet. 

Table 3. Targets with absorption systems showing marginal or significant variability. 


target 

4. 

date 

N 

cm“^ 

date 

N 

cm~^ 

AO0235-bl64 

Mg II 

(blend) 

marginal variability” 




0.85232 

12-1997 

(2.79 ± 0.46) X 10*2 

02-2007 

(2.82 ± 0.38) X 10*2 


0.85237'’ 

12-1997 

(5.58 ± 0.85) X 10*2 

02-2007 

(3.25 ± 0.57) X 10*2 


Mg II 


no variability 




0.85562 

12-1997 

(7.27 ± 0.36) X 10** 

02-2007 

(7.84 ± 0.28) X 10** 

PKS 1229-02 

Nai 


significant variability* 




-0.00013'' 

02-2002 

(1.98 ± 0.41) X 10** 

04-2011 

(4.6 ± 1.5) X 10*** 


-0.00010 

02-2002 

(1.81 ±0.20) X 10** 

04-2011 

(1.67 ± 0.16) X 10** 


Mg II 


no variability 




0.76836 

02-2002 

(7.15 ± 0.24) X 10** 

04-2011 

(7.57 ± 0.58) X 10** 


Mg II 


marginal variability” 




0.83032' 

02-2002 

(3.42 ± 0.07) X 10*2 

04-2011 

(2.87 ± 0.13) X 10*2 

FBQS J2340-0053 

Ci 

(blend) 

marginal variability* 




2.05470 

08-2006 

(1.33 ± 0.16) X 10*2 

10-2008 

(1.25 ± 0.07) X 10*2 


2.05473' 

08-2006 

(3.13 ± 0.49) X 10*2 

10-2008 

(1.52 ± 0.23) X 10*2 


H2,J=3 


marginal variability* 




2.05499 

08-2006 

(1.66 ± 0.10) X 10*2 

10-2008 

(1.43 ± 0.09) X 10*2 


2.05513' 

08-2006 

(6.19 ± 0.50) X lO*'* 

10-2008 

(4.39 ± 0.31) X lO*'* 


ajj.c,d ggg pjg 4 ^ 2, 5, 3 (left panel), 3 (right panel). 

Significant variability: significance level >3.5(T and AN/{N) > 25%. 

'' Marginal variability: significance level either in the range 2.0-3.5(T or >3.5(T but AN/{N) < 25%. 


der) expected to be regular, a third spectrum should confirm the 
trend observed on the two first ones. Such a work is currently 
in progress (see Appendix B2 in Dawson et al. 2013). Repeated 
high resolution observations of the brightest quasars involved 
would also be very useful to check the reality of time variations 
and to characterize them. 


4.3. High z neutral/molecular gas 


Ci and H 2 are known to be closely related in diffuse molec¬ 
ular clouds (Srianand et al. 2005). Indeed, the velocity com¬ 
ponents found for both species are generally similar, although 
some significant differences in their relative strength can be 
observed, as is the case for the Zabs = 2.0545 system toward 
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Fig. 4. Intervening Mgii absorption at Zabs = 0.852 toward 
AO 0235+164 in the December 1997 (HIRES, black line) 
and February 2007 (UVES, red line) spectra (lower panel: 
Mg 11/12796; upper panel: Mgii/12803). The blue vertical tick 
mark underlines a marginally variable blended component at 
2abs = 0.85237, although at a 2.3cr significance level only. 


FBQS J2340-0053 (Sect. 3.3 and Fig [3]). We successively dis¬ 
cuss these two tracers below. 

Ci lines are detected in three high z systems, one toward 
PKS 0458-020 and two toward FBQS J2340-0053, but useful 
constraints could be obtained only for the latter two at Zabs = 
1.3606 (no change in N) and 2.0454 (possible variation of about 
70% for an unresolved component). Unfortunately, since i) the 
transverse velocity of the quasar is unknown and ii) the LoS drift 
is determined solely by the quasar motion for the Zabs = 2.0545 
system because Zabs - Zem, it is not possible to compute the 
scale probed through the neutral gas. Further, the transverse ob¬ 
server’s velocity for this target being only 12 km s“', little drift 
due to the observer’s motion is expected for the Zabs = 1.3606 
absorber. We can then just get a typical value for the drifts by 
assuming a transverse velocity of about 300 km s“' for the in¬ 
tervening galaxies. This value together with a time interval of 
1.9 yrs (observer’s frame) between the two observations corre¬ 
sponds to a linear scale of about 50 and 40 au at Zabs = 1.3606 
and 2.0545 respectively. In our own Galaxy, the available infor¬ 
mation on structure in the C i distribution remains quite limited. 
Welty (2007) performed a detailed study of the time behavior 
observed for visible and UV transitions due to an intermediate 
velocity cloud toward the star HD 219188. Variations by a factor 
of about two are seen for Na i, C i and Ca ii column densities over 
scales of tens of au. Qualitatively, Nai and Ci display the same 
time behavior, indicating that the more extensive results obtained 
on the Nai structure (Crawford, 2003; Smoker et al. 2011) are 
probably valid also for C i. In this context, results obtained for 
V(Ci) at Zabs = 1.3606 and 2.0545 are certainly consistent with 



Fig. 6. S 11/1/11253,1259 absorption at Zabs = 2.05454 toward 
FBQS J2340-0053 in the november 2006 (HIRES, black line) 
and October 2008 (UVES, red line) spectra. 


the properties known for Galactic gas. Clearly, studies involving 
more systems are desirable to reach firm conclusions. 

Regarding H 2 , velocity components number 5 and 6 display 
an intermediate opacity which provides a good sensitivity for 
the 7 = 3 column densities (the other 7 = 3 components are 
either too weak or too saturated; similarly, 7 = 2 lines are too 
strong while 7 = 4 are too weak to be suitable). The detection of 
marginal variations in component number 6 is interesting in the 
context of questions related to the ’’warm H 2 ” observed along 
many Galactic LoS (Gry et al. 2002; Jensen et al. 2010). UV 
pumping is generally found to be insufficient to explain the ob¬ 
served excess of 7 > 2 H 2 with respect to the amount expected 
for a thermal distribution (Gry et al. 2002). Small structures (like 
shocks or turbulent vortices) where the temperature can locally 
reach high values have been proposed to account for both the 
production of species like CH^ and the relative amount of 7 > 2 
H 2 (Godard et al. 2009). Presently, there are no direct constraints 
on the size of these putative regions. In our Galaxy, the only re¬ 
peated measurements in the far-UV were performed toward the 
runaway O star, HD 34078 (Boisse et al. 2009). Stringent con¬ 
straints were derived from the analysis of damped line profiles 
associated with 7 = 0 lines, but no limit could be obtained for the 
7 = 2,3,4 levels because the corresponding lines are insensitive 
to small relative changes in N. It is noticeable that for the varying 
component toward FBQS J2340-0053, the N{J = 0,1,2,3) val¬ 
ues (1.39 xlO*^ 2.61 xlO*’, 4.71 xlO'^ and 5.28 xlO*'^ cm^^; 
for the 7 = 3 column density, we give the average of the 2006 
and 2008 values) are consistent with a unique excitation temper¬ 
ature, Tex = 120 K, unlike what is seen over most of Galactic 
LoS where an excitation temperature of about 300 K is required 
to fit N(J) values at 7 > 2 (Jensen et al. 2010). This indicates 
that the gas is dense enough to ensure thermalisation, in qualita¬ 
tive agreement with the small size suggested by the presence of 
time variations for this component. 
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4.4. Local Galactic gas 

We found a large variation (factor of four) in Galatic Ai(Nai) for 
one of the five velocity components toward PKS 1229-02. The 
constraint on the corresponding cloudlet size can be estimated 
from the peculiar motion of the Sun with respect to the LSR (see 
Sect. 2.2). For this direction, the Sun’s velocity is mostly trans¬ 
verse with V, (Sun)/LSR = 17.8 km s“\ implying a drift of about 
35 au (this is only an estimate since we cannot know the trans¬ 
verse velocity of the cloud). This is well within the range (typi¬ 
cally a few 10 au or larger) for which Na i line changes have been 
observed on other LoS (Crawford 2003). The other four compo¬ 
nents toward PKS 1229-02 remained stable over the same linear 
scale, indicating that the associated fragments are characterised 
by sizes larger than 100 au. For other targets providing useful 
constraints on Galactic gas structure, we find variation neither 
toward AO 0235-H164 nor toward FBQS J2340-0053. The trans¬ 
verse velocity of the Sun for these two targets is 15.8 km s“* and 
16.6 km s ' respectively, implying drifts of about 30 and 7 au 
in local material. In these two cases, Nat lines have a relatively 
large opacity implying a poor sensitivity to variations. 

Toward PKS 1741-038, the estimated drift over the 9.8 yr 
separating the two observations is about 17 au. The absence of 
variations in DIBs near 5780.5 and 5797.1 A can be related tore- 
suits obtained recently by Cordiner et al. (2013) for p Oph stars. 
These authors find variations of about 5-9% for the above DIBs 
over scales exceeding 300 au. Thus, given our limited sensitivity, 
the observed stability is not surprising. 

5. Summary and prospects for further studies 

As discussed above, the only clear and unambiguous variation 
has been detected for the Nat doublet from our own Galaxy to¬ 
ward PKS 1229-02, in good agreement with previous results on 
the au-scale structure of neutral Galactic gas. In the same kind of 
media but at higher redshifts (z ~ 2 toward FBQS J2340-0053), 
we also find tentative evidence for structure in neutral and 
molecular gas : C i with AN/{N) ^ 70% at the 3.0cr significance 
level and H 2 with AN/{N) ^ 35%, at the 3.1cr level. For low 
P absorption systems in blazar spectra, absorption profiles are, 
for the most part, very stable and only a marginal 3.7cr 17% 
variation has been seen in the strength of a narrow Mg ii compo¬ 
nent toward PKS 1229-02. On the other hand, no change was 
seen for moderately ionised gas at high p (i.e. intervening halo 
gas), which again is consistent with known constraints on typical 
cloud sizes in such media. 

The results discussed above illustrate well the potential and 
the limits of the method used in this paper to investigate the 
small-scale structure of absorbing gas in quasar spectra. Let us 
first outline the great value of high spectral resolution. Often, 
several transitions from the same species are detected and can 
be used to fit absorption line profiles, thus yielding well con¬ 
strained b and N values for the velocity components seen at each 
epoch. This redundancy provides a powerful objective way (at 
least for those components which are partially resolved) to as¬ 
sess the reality of time changes and rule out false variations. 
Indeed, it is very unlikely to find by chance two artefacts mim¬ 
icking true variations, i.e. occurring on two line profiles at the 
same position in velocity space and with a strength consistent 
with the / values of the transitions involved. If species display¬ 
ing more than two transitions (like Ci or I-I 2 ) can be used to 
increase the redundancy and the range covered by / values, the 
probability is even lower and better constraints can be derived 
from line fitting. Conversely, the Mg ii system at Zats - 0.85236 


toward AO 0235 H-164 illustrates well that with only two tran¬ 
sitions differing in their / values by a factor no larger than two, 
one may be unable to reach firm conclusions when variations are 
moderate. Furthermore, high resolution allows us to quantify the 
changes in meaningful astrophysical parameters, N and b values, 
or to set upper bounds on their variations. Limitations appear 
however for some line profiles (e.g. the Ci component toward 
FBQS J2340-0053 for which variations are suspected), when 
the decomposition into discrete components required by line fit¬ 
ting is poorly defined or not unique. Such difficulties mean that 
the spectral resolution used is not high enough to analyse prop¬ 
erly the velocity distribution, implying that line opacities as well 
as N and b values cannot be determined reliably. One should 
therefore prefer absorption systems with simple velocity struc¬ 
ture for repeated observations. 

In our study, we considered several topics for which useful 
constraints can be obtained from an investigation of time varia¬ 
tions. We now discuss prospects for two of them. We mentioned 
above that generally, we have a poor knowledge of the scales 
probed by repeat observations because the transverse motion of 
the target is not measurable. The situation is more favorable for 
low redshift absorbers because the drift of the LoS through the 
intervening gas is determined mainly by the observer’s motion, 
which is inferred directly from observations of the CMB dipole. 
For instance, in order to study whether the small-scale structure 
seen in the Galactic neutral or weakly ionised gas is ubiqui¬ 
tous in external galaxies, one could reobserve quasars belong¬ 
ing to the so-called quasar-galaxy pairs (see e.g. Bowen et al. 
1991) and search for temporal changes in the Can or Nat lines 
over more than 20 years. Of course, targets located in directions 
which are more or less normal to our large-scale peculiar veloc¬ 
ity should be preferred (in this favorable configuration, the scales 
probed over 20 years are of the order of 1500 au). Regarding 
high redshift absorbers, for which only a rough estimate of the 
scales probed can be obtained, one needs to observe a statisti¬ 
cally significant sample of targets to derive meaningful results 
concerning the presence and strength of small-scale structure. 
Some quasars are found in pairs (Hennawi et al. 2010) and occa¬ 
sionally, absorbers appear to belong to galaxy clusters (Whiting 
et al. 2006) implying increased peculiar velocities. Such LoS 
will, in average, be associated with larger drifts, and should be 
targetted preferentially. 

A second topic for which repeated observations of quasars 
at high spectral resolution are especially timely involves the 
properties and origin of the warm (7 > 2) H 2 found in diffuse 
molecular gas in our own Galaxy, as well as in high redshift 
absorbers. During the past ten years, the number of high z H 2 
systems has increased very significantly and now exceeds 20 
(Albornoz-Vasquez et al. 2014). Some of them display narrow 
velocity components which are well suited to a detailed study 
(see e.g. Noterdaeme et al. 2007). As discussed in Sect. 3 for 
the system at Zabs = 2.0545 toward FBQS J2340-0053, multi¬ 
epoch observations can provide key information on the typical 
size of regions where diffuse molecular gas is heated to temper¬ 
atures above 200 K. Such high z warm H 2 is seen at epochs when 
star formation was much larger than today (Madau & Dickinson 
2014) and one can therefore anticipate that the diffuse molecular 
medium should be strongly impacted by a very high dissipation 
rate of the mechanical energy carried by shock waves and turbu¬ 
lence. 

To conclude, we wish to stress that the data used here were 
acquired for purposes other than variability searches and thus, 
are certainly not optimal for our study. First, the two spectra 
available for each target often differ significantly in their S/N 
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ratio and one is then limited in the comparison by the lower 
quality spectrum. The time interval also varies a lot from target 
to target (unfortunately, the lag for the quasar displaying the 
richest spectrum - FBQS J2340-0053 - is of only two years). 
A straightforward strategy to get better constraints or address 
some questions discussed above with more appropriate data 
would consist in i) selecting targets displaying the required 
types of absorption lines and with one already available good 
high resolution spectrum taken at least 10 years ago, and ii) 
reobserving them at a similar S/N and resolution. Obviously, 
the existence of databases where high resolution quasar spectra 
are archived properly are extremely useful in such an approach. 
Finally, we note that improvements in the accuracy of wave¬ 
length calibrations accomplished recently in the context of 
the ESO Large Program designed to constrain the variation 
of fundamental constants (Molaro et al. 2013) or forseen for 
instruments that will become available in the near future (e.g. 
the ESPRESSO spectrograph to be mounted on the ESO/VLT 
telescope) should indeed help to investigate in more detail the 
velocity field within the absorbers. 
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